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Abstract 

We study discrete vortices in coupled discrete nonlinear Schrodinger equations. We 
focus on the vortex cross configuration that has been experimentally observed in pho- 
torefractive crystals. Stability of the single-component vortex cross in the anti-continuum 
limit of small coupling between lattice nodes is proved. In the vector case, we consider two 
coupled configurations of vortex crosses, namely the charge-one vortex in one component 
coupled in the other component to either the charge-one vortex (forming a double-charge 
vortex) or the charge-negative-one vortex (forming a, so-called, hidden-charge vortex). We 
show that both vortex configurations are stable in the anti-continuum limit if the param- 
eter for the inter-component coupling is small and both of them are unstable when the 
coupling parameter is large. In the marginal case of the discrete two-dimensional Man- 
akov system, the double-charge vortex is stable while the hidden-charge vortex is linearly 
unstable. Analytical predictions are corroborated with numerical observations that show 
good agreement near the anti-continuum limit but gradually deviate for larger couplings 
between the lattice nodes. 



1 Introduction 

In the past few years, the developments in the nonlinear optics of photorefractive materials 
m and of Bose-Einstein condensates in optical lattices 121 01 ^ have stimulated an enormous 
amount of theoretical, numerical and experimental activity in the area of discrete nonlinear 
Hamiltonian systems. A particular focus in this effort has been drawn to the prototypical 
lattice model of the discrete nonlinear Schrodinger (DNLS) equation [H]. The latter, either as 
a tight binding limit , or as a generic discrete nonlinear envelope wave equation plays a 
key role in unveiling the relevant dynamics within the appropriate length and time scales. 

One of the principal directions of interest in these lattice systems consists of the effort to 
analyze the main features of their localized solutions. In the particular case of two spatial 
dimensions, such structures can be regular discrete solitons 8 or discrete vortices (i.e., struc- 
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tures that have topological charge over a discrete contour) The study of these types of 
coherent structures has made substantial leaps of progress in the past two years with the nu- 
merical and experimental observation of regular discrete solitons ^l^J, dipole solitons ;1.2j . 
soliton-trains ^Hj , soliton- necklaces and vector solitons ^3] in photorefractive crystals and 
experimental discovery of robust discrete vortex states ^lEl; based on earlier theoretical 
predictions p ilTIH Onj. 

On the other hand, the recent years were marked by the experimental developments in soft 
condensed-matter physics of Bose-Einstein Condensates (BECs). Among the important recent 
observations one can single out the experimental illustration of the dark |21l 1221 123j , bright 
|241 125j and gap I^H] solitons in quasi-one dimensional BECs. The experimental capabilities 
seem to be on the verge of producing similar structures in a two-dimensional context |27j . 

In both of the above contexts (nonlinear optics and atomic physics), multi-component systems 
were recently studied due to their relevance to applications. In particular, the first observa- 
tions of discrete vector solitons in nonlinear waveguide arrays were reported in |281 129j . while 
numerous experiments with BECs were directed towards studies of mixtures of different spin 
states of s^Rb (SOI EI] or '^^Na jS2] and even ones of different atomic species such as ^^K-^'^Rb 
|33j and ^Li~^^^Cs While the above BEC experiments did not include the presence of an 
optical lattice, the addition of an external optical potential could be manufactured within the 
present experimental capabilities 0. 

It is the purpose of the present work to address these recent features of the physical experi- 
ments, namely discrete systems with multiple components. In particular, we aim at addressing 
the fundamental issue of how localized excitations are affected by the presence of two com- 
ponents which are coupled (nonlinearly) to each other. While our results will be presented 
for the specific example of two coupled DNLS equations with cubic nonlinear ities, we believe 
that similar features persist in a variety of other models. We should note here that rather 
few studies have focused on the two-dimensional vector generalization of the DNLS equation 
|351 1361 157] . Among others, we mention the work j38j which was motivated by the experimen- 
tal system of the nonlinear waveguide arrays proposed in j^H]. To the best of our knowledge, 
these earlier studies did not address vortices in coupled discrete systems. 

For vortices in coupled systems, a number of interesting questions emerges concerning the 
stability of particular vortex configurations (e.g. the so-called vortex cross ^JE]) including 
the case of equal charges in both components and the case of opposite charges between the 
two components. The former state has a double vortex charge, while the latter has a hidden 
vortex charge. It has been shown for the continuous NLS equation with cubic-quintic [22] and 
saturable jlU] nonlinearities that these two states have different stability windows. 

In the present setting, we examine the stability of such vortex structures in the discrete 
case both analytically and numerically. We use the method of Lyapunov-Schmidt reductions 
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developed earlier in [nj. This method allows for direct analytical calculations of eigenvalues 
of the linear stability problem as functions of the system parameters (such as the coupling 
between adjacent lattice sites and the coupling between the two components). 

Our presentation is structured as follows. In section 2, we introduce the setup and the vor- 
tex cross configurations. In section 3, we study the stability of such configurations in the 
one-component model. In section 4, we generalize the vortex cross configuration to the two- 
component case and compare our results with numerical computations of the parameter con- 
tinuations. In section 5, we deal with a special Manakov case of the system of two DNLS 
equations. Finally, in section 6, we summarize our findings. Appendix A presents technical 
details for the case of the single-component vortex cross. 

2 Setup 

We write the coupled system of discrete nonlinear Schrodinger (DNLS) equations in the form: 

^Ufi^m ~l~ £ (^n+l,m ~l~ ^n— l,m ~l~ '^n.m+l ~l~ Un^m—l} ~l~ (|''^n,m| ~l~ /?|f^ra,m| ^Un^m — 0) (^'1) 
^^n,m ~l" £ i'^n+l,m ~l~ ^n— l,m ~l~ ^n,m+l ~l~ f^ra,m— l) ~l~ (/5|^n,m| ~l~ |f^n,m| )^n,m — 0) (^•^) 

where /3 is a non-negative parameter for the coupling between the two components (n, v) and 
e is a small non-negative parameter for the coupling between adjacent lattice sites. Localized 
modes of the coupled system ()2.1() - H2.2() take the form: 

where a; is a parameter of time-periodic solutions and {4>n,rm'^n,m) satisfy the system of non- 
linear difference equations: 

(1 - \(t>n,m\'^ - P\lpn,m\'^)4>n,m = £ {4>n+l,m + 4>n-l,m + <t>n,m+l + <Pn,m-l) , (2.4) 
{lV - (3\(t)n,m\^ - \ll^n,m\^)'4'n,m = £ (V'n+l.m + ^n-l,m + V'n.m+l + V'n,m-l) • (2.5) 



We are interested in a particular vortex solution, called the vortex cross. An example of this 
solution is obtained numerically for /3 = |, w = 1 and e = 0.1 and it is shown on Figure ^ 
Let us consider the diagonal square discrete contour on the grid (n,m) G 1?: 

5(0) = {(-1,0); (0, -1); (1, 0); (0, 1)} C I?, (2.6) 

enumerated in the same order by j = 1, 2, 3, 4. We shall assume that the vortex cross of Figure 
n bifurcates from the limiting solution at the anti-continuum limit e = 0: 

(0)^^ae*^^ (n,m)G5W ,(o)^^^e-^ (n,m)G5W 

1 0, (n,m)^5(0) "^"'"^ I 0, (n,m)^5(0) ^ ' 
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Figure 1: The contour plots show the amphtude and phase (left and right panels respectively) 
of the two components (top and bottom respectively) for a (1, 1) (left four subplots) and a 
(1, —1) (right four subplots) vortex configuration, in the case of /3 = 2/3, lj = 1, and e = 0.1. 



where the set of phase parameters {9j, i^j}j=i is yet to be determined, while the set of amplitude 
parameters (a, b) is determined from solutions of the system: 



(2.i 



(2S 



+ = 1, (3a^ + = uj. 
When /? / 1, there exists a unique solution of the system ()2.8|) : 

2 _ 1- Puj 2 _ P 

The solution is meaningful only if > and 6^ > 0, which define the domain of existence: 

min(/3, (3^^) < to < max(/3, jS'^). (2.10) 

When /3 = 1, the domain of existence shrinks into the line w = 1 and the solution of the 
system (|2.8j) forms a one-parameter family: 



a = cos 6, b = sin 6, 6 £ [0, 2tt]. 



(2.11) 



The vortex cross, if it exists, is defined by the phase configurations along the discrete contour 
5(0): 

7r(i - 1) Tr{j - 1) 
^ ' - ± n ' 



j = 1,2,3,4. 



(2.12) 



The upper sign corresponds to the (1,1) coupled state called the double-charge vortex, while the 
lower sign corresponds to the (1,-1) coupled state called the hidden- charge vortex. Persistence 
and stability of the vortex configurations ()2.7|) . ()2.9j) . and 1)2. 12() are addressed separately in 
the cases /3 = 0, < /3 < 1, /? = 1, and /? > 1. 
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3 Scalar vortex cross 



We apply the method of Lyapunov-Schmidt (LS) reductions developed in W\ to the scalar 
nonlinear difference equation: 

(1 - \^n,m\^)'^n,m = £ (^n+l,m + ^n-l,m + '^n,m+l + '^n,m-l) , (3.1) 

This scalar equation corresponds to the reduction tpn,m = 0, V(n, m) £ H? of the system 
()2.4I) - H2.5() . Local existence of a single-component vortex cross in the scalar problem 1)3. 1|) is 
proved in Appendix A for small values of e (on the basis of Proposition 2.9 in ^). This result 
is formulated as follows. 

Proposition 3.1 There exists a unique (up to the gauge invariance) continuation in e of the 
limiting solution at e = 0: 

$(0) =h'''^ (n,m)ES(0) 
"•"^ \ 0, (n,m)^5W ^ ^ ^ 

where S^'^^ is given by i2.b]) and the values of 6j are given by \2.1^) . The family of vortex 
solutions $n,m(e)) {n,m) £l? is a smooth (real analytic) function of e. 

To address spectral stability of the vortex cross in the time-evolution of the single-component 
DNLS equation, we consider the linearization problem with the explicit formula 

'Un,m(^) — 6 ^n,m ~\~ an,m(^ ~l~ bn^m^^ 

and derive the linear eigenvalue problem from the DNLS equation, 

(1 2|<I'n,m| )an^m ^ n,mbn,m ^ian+l,m, ~\~ On— l,m ~l~ On^rra+l ~l~ O-n^m- 1) — i^an^m (3. 3) 
= 2 



m^n,m ^{bn+l,m ~l~ ^n— l,m + "n,m+l iXbn^m, (3.4) 



where A is an eigenvalue and {an,m-, bn,m) are components of an eigenvector. Symbolically, we 
write the linear eigenvalue problem as 

'H{e)ip = iXaif, (3.5) 

where 'H(e) is the linearized Jacobian matrix for the system ()3.1() . o" is a diagonal matrix of 
(1,-1), and (p is an eigenvector consisting of {an,m,bn,m)- The linear eigenvalue problem for 
the limiting solution $n,m = ^n}n at e = has a set of double zero eigenvalues with the 
eigenvectors ej and generalized eigenvectors e^, such that TC^^^ej = and TC^'^^ej = 2iaej, 
where Ti^'^^ = 7i{0). The index j enumerates the set S^^^ and the eigenvectors and ej have 
non-zero components only at the corresponding nodes of the set S^^^ , 

(3.6) 









( e'^^ \ 


ej = i 
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The kernel of 'H(e) for e 7^ includes at least one eigenfunction 

(n,m)GZ^ (3.7) 

which follows from the gauge invariance of the DNLS equation with respect to rotation of 
the complex phase in ^n,m, {n,m) E 1? . It is easy to show that a generalized kernel for zero 
eigenvalue is non-empty as it includes a solution of the inhomogeneous equation TC{e)(p = 2ia(p 
exists, where (f is given by (|3.7jl . 

Using the perturbation series expansion for $n,m(e), we define the expansion 'H(e) = TC^^^ + 
eTC^^^ + e^Ti^"^^ + O(e^). By Lemma 4.1 in 9 , computations of Appendix A determine the 
splitting of zero eigenvalues of 'H(e) as e / 0. The splitting of zero eigenvalues of a7i{e) is 
formulated and proved as follows. 

Proposition 3.2 Let ^n,m{^), {n,m) G be a family of vortex solutions defined by Propo- 
sition 15'. il The linearized problem ici. has zero eigenvalue of algebraic multiplicity 
two and geometric multiplicity one and three small pairs of purely imaginary eigenvalues of 
negative Krein signatures^ with the asymptotic approximations, 

Ai,2, A3,4 = ±2ie + 0(e2), Ag.e = ±^ie^ + ©(e^). 

The rest of the spectrum is bounded away the origin as e — > and it is located on the imaginary 
axis of A. 

Proof. We supplement the general proof of Lemma 4.2 in |9| with the explicit perturbation 
series expansions for small eigenvalues of the linear eigenvalue problem (|3.5|) : 

<^ = (^(0) + e(^« + eV^'^ + O(e='), A = eAi + e^Aa + 0(e3), (3.8) 

where 

^ — Z^h^]-' ^ — 2 / > J J ^ tinhorn' 

and the solution <J{Xo^ = -^("^"^^(^V^"'' = -W^^V^"-* is uniquely defined on the set 5^^), 
where S^^^ is the set of adjacent nodes to the set of S^'^h At the second-order perturbation 
theory, the problem is written in the form, 

7^(0)^^(2) + ^(1)^(1) + ^(2)^(0) ^ -^^^^H) ^ ix^a^p^'^l (3.9) 
Projecting the problem to the kernel of TL^^\ we find the reduced eigenvalue problem: 

A^2C=^A?c, (3.10) 

simple eigenvalue of the linear eigenvalue problem l|3.5^ is said to have the negative Krein signature if 
the quadratic form for the associated eigenvector {T-l{e)(fi, ip) is negative. 
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° 0.1 0.2 e 0.3 0.4 

Figure 2: Eigenvalues of the scalar vortex cross versus e. The top panel shows the imaginary 
part of the relevant eigenvalues, while the bottom panel shows the real part. The solid lines 
display the numerical results, while the dashed ones correspond to the asymptotic approxima- 
tions. 

where c = (ci, 02,0^,04)'^ and M-2 is computed in Appendix A. Therefore, two negative eigen- 
values 71^2 of the Jacobian matrix e^A^2 generate two pairs of imaginary eigenvalues of negative 
Krein signatures in the linear eigenvalue problem by virtue of the relation A = The 
same computation is then extended up to the fourth order, where it is found that the nega- 
tive eigenvalue 73 of the extended matrix e^7W2 + e^A^4 determines the third pair of purely 
imaginary eigenvalues by virtue of the same relation A = ■ 

We note that the count of eigenvalues of negative Krein signatures corresponds to the closure 
theorem for negative index of 'H(e) (see for details). There are four negative eigenvalues 
of 7i^^^ for the limiting solution H3.2|) and three more small negative eigenvalues occur for 
e 7^ 0. The total number of negative eigenvalues is reduced by the gauge symmetry constraint, 
such that six negative eigenvalues in a constrained subspace match three pairs of imaginary 
eigenvalues with negative Krein signature. 

The asymptotic approximations of eigenvalues A are plotted on Figured by dashed lines. The 
numerical computations of the same eigenvalues (up to the prescribed numerical accuracy) 
versus e are shown by solid lines. All three pairs of purely imaginary eigenvalues bifurcate 
into complex domain when they collide to other eigenvalues of stability problem (e.g. with 
eigenvalues of positive Krein signatures or with the spectral band). The first collision is 
numerically detected to occur at e ~ 0.395. 
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4 Vector vortex crosses for < < 1 and (3 > 1 



In order to consider the coupled vortex configurations in the non-degenerate case /? 7^ 1, 
we extend computations of Appendix A to the solution of the coupled nonlinear difference 
equations H2.4() - ()2.5() . We report here computations for two related problems: (i) bifurcations 
of small eigenvalues of the linearized Jacobian matrix near the zero eigenvalue and (ii) bi- 
furcations of small eigenvalues of the linearized stability problem near the origin. Because 
of the computational complexity of the analytical approximations, we shall complement the 
analytical results of the second-order Lyapunov-Schmidt (LS) reductions with the symbolic 
computational results of the fourth-order LS reductions. 

Similarly to the scalar case, the linearized stability problem for the two-component system 
takes the matrix-vector form: 

n{e)(p = iXaif, (4.1) 

where 'H(e) is the linearized Jacobian matrix for the system (|2.4() - (|2.5() . a is a diagonal matrix 
of (1, —1, 1, —1), and (p is an eigenvector consisting of four elements of the perturbation vector 
at each node {n,m) G Z^. The diagonal block of the matrix 'H(e) at each node {n,Tn) £ 1?' 
takes the form: 

/ 1 - 2|(/)„^TOp - /3|'(/'„,mP m -P4'n,m^n,m -P4'n,m^n,m \ 

-P(t>n,m'lpn,m - l^4>n,mi^n,m W - /3| 0„,m ^ - 2 j-^n^m P -^n,m 

The non-diagonal blocks of 'H(e) comes from the difference operators in the right-hand-side of 
the system ^I ^ -^I ^ . 



4.1 Bifurcations of zero eigenvalues of the linearized Jacobian matrix 

We extend the perturbation series expansions ()A.1|) to the two-component case, 

oo oo 

^n,M = Yl ^"'^S^' ^n,m{e) = J] e^S^L, (4.2) 



fc=0 



fc=0 



where the zero-order solution in the anti-continuum limit is given by (|2.7|) and parameters 
(a, 6) are given in 1)2. 9() . The first-order corrections are found from the uncoupled system of 
equations, similarly to the scalar case: 



' 0, (n,m) G 5(0) 



n, m 



g5(i) 



v^(i) 



^ 0, (n,m) ^ 5(0) U 5(1) 



' 0, (n,m) G 5(0) 
oj-^b^l^K'''^ (n,m)G5(i) (4.3) 
0, (n,m) ^ 5(o)u5(i) 
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where is defined in (|A.3() . The second-order corrections are found in the form: 



r'n.m 



' sfe'^^, (n,m) G 5^ 
0, (n,m)G5(i) 

0, (n,m) ^ 5(0) U5«U 5(2) 



(2) 
n,m 



0, (n,m)E5(i) 



0, (n,m) ^ 5(0) U5(i) U5(2) 

(4.4) 

where X^P^ is defined in HA.5|) . The real parameters satisfy an inhomogeneous 

system 



2{asf +Pbrf) 



2{(3asf +br'f>) 



{2)^ 



4 + 2 cos(6lj+i - Oj) + 2 cos(6lj_i - 6lj) + cos(ej+2 - Oj), (4.5) 
4 + 2 cos(z^j+i — Vj) + 2 cos(z/'j_i — i^j) + cos(z/'j4.2 — Vj)- (4.6) 



When /3 / 1, the inhomogeneous system ()4.5() - (|4.6() has a unique solution. Second-order 
corrections to the bifurcation equations are uncoupled and have the form: 



a 



(2) 

i 



2sin(6'j - 6j+i) + 2sin(6'j - ^^i) + sin(6'j - 6*^+2), 
2sin(z^j - Vj+i) + 2sin(z^j - Vj-i) + sin(i/j - i'j+2)-, 



(2) (2) 

where a suitable normalization of g)j and is made. As a result, the Jacobian matrix 
computed from derivatives of {g^p , ^^^'^ (^«' block-diagonal as diag(A42, ^^2)1 where 
A^2 is defined in Appendix A. By Lemma 4.1 in j^j, non-zero eigenvalues of diag(7W2, -^2) 
determine small eigenvalues of the linearized Jacobian matrix 'H(e), 

71,2,3,4 = -2e2 + 0(6^). 

Two zero eigenvalues of diag(A^2) -^^2) split into two non-zero eigenvalues in the fourth-order 
LS reductions, while two other zero eigenvalues of diag(A^2i -^^2) persist beyond all orders due 
to the gauge invariance of each component in the the coupled DNLS equations H2.1() - (|2.2I) . 
Indeed, the kernel of TC^e) for e 7^ includes at least two eigenfunctions: 





/ 




/ 







< 




(pn,m 
























v 


/ 


v 




) 



(n, m) € 



(4.7) 



In order to compute the small non-zero eigenvalues of the linearized Jacobian matrix 7i{e), 
we use the symbolic computation package based on Wolfram's Mathematical. The projection 

■^The software programs are available online at|littp://dmpeli.math. mcmaster.ca/Software/LSreductions.litmlj 
For more information on the symbolic mathematics package in which these programs were implemented, see 
|http://www.wolfram.com'] 
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to the eigenspace of diag(A^2, -^^2) spanned by eigenvectors (p2,04)^ and (04,p2)"^, where 
P2 = (—1, 1, —1, 1) and O4 = (0, 0, 0, 0), leads to the reduced eigenvalue problem (for m = 1), 

-8 



(ai ± Pa2) 



l + P 

{±Pai + 02) 



7"! 



702, 



where (ai, 02) are coordinates of the projections, 7 = linie^o £""^7; and the upper/lower signs 
refer to the two coupled vortices (l,ibl). It is clear that the eigenvalues of the reduced 
eigenvalue problem are the same for either sign and they define two small eigenvalues of the 
linearized Jacobian matrix 'H(e) (for uj = 1): 



75 



-8e^ + 0(e''), 



76 



(l + (9) 



4.2 Bifurcations of zero eigenvalues of the linearized stability problem 

We consider the eigenvalue problem (|4.H) in the limit of small e. Let TC^e) = TC^^^ + eTC^^^ + 
€^^(2) _^ 0{e^). The set of eigenvectors of H^^'^ej = and W(°)fj = takes the form: 



( 



\ 



e 3 



\ 








/ 



V 



(4.8) 



The corresponding set of generalized eigenvectors of TC^'^^ej = 2iaej and Ti.^'^^ij = 2iafj takes 
the form: 

/ A+e'^3 \ ( A.e'^3 \ 



A+e-^^3 



where 



A, 



B, 



A^ 



-P 



A^e~'^3 
B^e-^''^ 

B 



(4.9) 



a2(l-/32)' — ab{l-P^)' 62(1-/32)' 

Bifurcations of zero eigenvalues of the linear eigenvalue problem (|4.H) can be computed with 
the extended perturbation series expansions 1)4. 2() for (/>„,m(e) and ipn,m{^) and extended per- 
turbation series ()3.8|) for (p and A, where 



inhom ' 



i=i 



i=i 



j=i 
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and cp^^^^^ = -W(°)"iW(iV^°^ = -W(^V^°^ is uniquely defined on the set S^^l At the 
second-order perturbation theory, we have the same problem H3.9() . from which we derive the 
reduced eigenvalue problem: 

M2C = ^Xl{A+c + A.d) (4.10) 
M2d = ^A? (B+c + 5_d) , (4.11) 

where c = (ci, C2, C3, 04)-^, d = (di, ^2,^3, ^4)"^, and Ai2 is the same as in the scalar case. Let 
7i = ^Xf. The reduced eigenvalue problem (|4.1U|) ~ (|4.11|) has four zero roots for 71 and two 
double-degenerate non-zero roots for 71, given from the quadratic equation: 

(71 + 2a2)(7i + 26^) = iaH^p\ (4.12) 

If (J = 1, such that = 6^ = j^, then the two non-zero roots for 71 are found explicitly, 

2(lT|/3|) 

By using the relation Ai = ziz^/2y[, we have just proved that the linear eigenvalue problem 
()4.1|) in the case oj = 1 and < /5 < 1 has four small pairs of purely imaginary eigenvalues 
with asymptotic approximations: 



Ai,2, A3,4 = ±2ie + 0(e2), Ag.e, X7,8 = ±2ie^ ^--^ + 0{e^). 

Two pairs of eigenvalues X^^ and A7^8 become pairs of real eigenvalues in the case (3 > 1. 
Two pairs of zero eigenvalues of the reduced eigenvalue problem (|4.10|) - (|4.11j) split at the 
fourth-order LS reductions as pairs of non-zero eigenvalues Ag^io and Aii^i2. Two other pairs 
of zero eigenvalues persist beyond all orders for e 7^ since the geometric kernel includes 
two explicit solutions (UZJ) ^-ud there exists a two-parameter solution of the inhomogeneous 
equation 7i{e)ip = 2ia(p, where (p is given by ()4.7() . In order to find the small non-zero 
pairs of eigenvalues, we apply again the symbolic computation package based on Wolfram's 
Mathematica. The projection to the eigenspace of diag(Al2 5 -^^2) spanned by eigenvectors 
(p2,04)^ and (04,p2)"^ for Ai = leads to the reduced eigenvalue problem (for uj = 1), 

-8 1 2 

(ai±/3a2) = xt; ^A2(ai - /3a2) 



1 + ' 2(1-/3) 

-8 1 2 

(±/Jai-|-a2) = -^Tf. ^A2(-/3ai + 02), 



1 + /?^ ' ^ ^' 2(1-/3) 

where (ai,a2) are coordinates of the projections and the upper/lower signs refer to the two 
coupled vortices (1, ±1). The eigenvalues of the reduced eigenvalue problem differs between the 
double-charge vortex (1,1) and the hidden-charge vortex (1,-1). For the double-charge vortex. 
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the two pairs of small eigenvalues of the linearized stability problem are purely imaginary for 

any /3: 

1-/? 
l + P 

For the hidden-charge vortex, the two pairs of small eigenvalues of the linearized stability 
problem are purely imaginary for < /3 < 1 and real for /? > 1: 



(1,1): A9,io = ±4ie2 + O(e3), Aii,i2 = ±4z 



+ 0(e3). 



(1, -1) : A9,io, An,i2 = ±4i^ l-p|e2 + 0{e'). 

We can specify precisely how many purely imaginary eigenvalues of the linearized stability 
problem H4.ll) have negative Krein signature. When < (3 < 1, there are eight negative 
eigenvalues of 7i^^^ for the limiting solution 1)2. 7() and six more small negative eigenvalues 
occur for e 7^ 0. The total number of negative eigenvalues is reduced by two gauge symmetry 
constraints, such that twelve negative eigenvalues in a constrained subspace match six pairs 
of imaginary eigenvalues with negative Krein signature. When /? > 1, there are four negative 
eigenvalues of TC^'^^ for the limiting solution (|2.7|) and five more small negative eigenvalues 
occur for e 7^ 0. The total number of negative eigenvalues is reduced by one^, such that eight 
negative eigenvalues in a constrained subspace match two real eigenvalues and three pairs of 
imaginary eigenvalues with negative Krein signature for the double-charge vortex and four 
real eigenvalues and two pairs of imaginary eigenvalues with negative Krein signature for the 
hidden-charge vortex. Therefore, the last pair of purely imaginary eigenvalues Aii^i2 for the 
double-charge vortex has positive Krein signature for /? > 1. 

We obtain numerically small eigenvalues A for small values of e and lo = 1. The results are 
shown on Figure El for /3 = I and on Figure |1] for (3 = 2. The left plot corresponds to the 
vortex pair (1,1), while the right plot corresponds to the vortex pair (1,-1). We note that 
the degeneracy of the pairs Ai^2 = A3^4 and X^fi = Ar^s is preserved for the case (1, —1), such 
that each bolded curve is double. The degeneracy of these eigenvalues is broken for the case 
(1, 1) and it is also broken for the pair Ag^io 7^ An. 12 for the case (1, —1). 

In the case of /5 = |, shown in Fig. |21 all six pairs of neutrally stable eigenvalues bifurcate 
to the complex plane for larger values of e due to the Hamiltonian-Hopf (HH) bifurcation. 
The first HH bifurcation happens earlier for the case (1, 1) at e ~ 0.395, due to the broken 
degeneracy between the two pairs of eigenvalues Ai^2 and A3^4. For the case (1, —1), the first 
HH bifurcation occurs at e ~ 0.495, i.e. the the hidden-charge vortex has a larger stability 
window for < /? < 1 (a similar observation is reported for continuous systems in |39[ I40j ). 

In the case of /3 = 2, shown in Fig. ^ both cases (1, 1) and (1,-1) are always unstable due 

to the pairs of eigenvalues As^g and A7_8- There are also additional observations. In the case 

■^When /3 is increased from /3 < 1 to /3 > 1, the Hessian matrix related to two gauge symmetry constraints 
loses one positive eigenvalue that passes through zero at /3 = 1 to the negative eigenvalue for /3 > 1 |41| . 
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E E 

Figure 3: Eigenvalues of the vector vortex cross with uj = 1 and /3 = | versus e. Left: 
(1,1). Right: (1,-1). The sohd hnes show the numerical results, while the dashed lines 
show the asymptotic approximations. Bolded curves correspond to double eigenvalues (that 
remain indistinguishable within the parametric window examined herein). A good agreement 
is observed for e < 0.1. 



(1, 1), the pairs of double real eigenvalues in the second-order LS reductions X^^ and Xj^s split 
as a quartet of complex eigenvalues, similarly to our computations in jOj. Real and imaginary 
parts of the quartet of complex eigenvalues are shown on Fig. 4 (left panel) by bolded curves. 
Only three HH bifurcations out of four pairs of purely imaginary eigenvalues occur for larger 
values of e. In the case (1, —1), two more pairs of real eigenvalues occur such that the hidden- 
charge vortex is more unstable compared to the double-charge vortex for /? > 1. Only two HH 
bifurcations occur for large values of e. 

5 Vector vortex cross for (3=1 

In the case /? = 1, the existence domain of the coupled vortex configurations shrinks to the line 
UJ = 1. The zero-order solution in the anti-continuum limit is given by (|2.7() . where parameters 
(a, b) are given by (|2.11|) . The second-order solution of the linear inhomogeneous system (|4.5() - 
(|4.6|) with a singular matrix exists provided that the values of 6j and z/j are defined by (|2.12l) . 

(2) (2) 

The arbitrary parameter in the second-order solution Sj and renormalizes the arbitrary 
parameter 6 in the representation (|2.11|) . 

When (3 = UJ = 1, the existence problem (|2.4() - (|2.5() is symmetric with respect to components 
i4'n,m,''Pn,m) such that the system (|2.4() - ()2.5() can be reduced to the scalar difference equation 
(|3.1|) with the two independent transformations: 

(1,1) : 4>n,m = cos (5 ^n,m, i^n,m = sin (5 (^n,m, 
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Figure 4: Eigenvalues of the vector vortex cross with uj = 1 and /? = 2 versus e. Left: 
(1, 1). Right: (1, —1). The sohd hnes show the numerical results, while the dashed lines show 
the asymptotic approximations. Bolded curves on the left panel correspond to the real and 
imaginary parts of complex eigenvalues, while bolded curves on the right panel correspond to 
double eigenvalues. 



(1,-1) : (j)n,m = cos 5 ^n,m, ^pn,m = s'm 5 ^n,m- 

The existence result for the scalar vortex cross is formulated in Proposition We will need 
the following non-degeneracy condition for the scalar vortex cross: 

E ^Irn ] ■ (5.1) 

It is clear from the limiting solution 1)2. 7() that the constraint 1)5.11) is satisfied for small e. The 
stability problem (|4.1|) is different between the cases (1, 1) and (1,-1). 




5.1 Eigenvalues of the (1,1) vortex cross 

In this case, the stability problem ()4.1)) is block-diagonalized under the following transforma- 
tion of the four components of the vector (p on the lattice node {n,m) G Z^: 



(1,1) 



/ 






/ 


cos 5 





sin 5 











bn,m 









cos 5 





sin 5 








C+ 






— sin (5 





cos 5 









V 




) 


V 





— sin (5 





cos 5 


) 


\V4 ) 



The components {an,m-,bn,m) satisfy the linear eigenvalue problem for scalar vortices 1)3. 3|1 - 
()3.4|) . The components (c^^rw ^n,m) satisfy two uncoupled self-adjoint eigenvalue problems: 
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Using the result of Proposition 13 . 21 and equivalent computations for the uncoupled self-adjoint 
problems ()5.2|) . we prove the following result. 



Proposition 5.1 Let $n,m(e); {n,m) G be a family of vortex solutions defined by Propo- 
sition W^ The linearized problem i4-l\ ) case /3 = 1 for the (1, 1) vortex cross has zero 

eigenvalue of algebraic multiplicity six and geometric multiplicity five and five small pairs of 
purely imaginary eigenvalues given asymptotically by 

Ai,2, A3,4 = ±2ie+0(e2), A5,6 = ±2ie^+0{^), Ay.s = iGie^+oCe^), Ag.io = ±^i€^+0{€^). 

The rest of the spectrum is bounded away the origin as e ^ and it is located on the imaginary 
axis of A. 

Proof. It remains to study bifurcations of zero eigenvalues in the self-adjoint problem ()5.2() 
as e 7^ 0. Let us define the perturbation series for the problem ()5.2I) : 

C± = c(0) + ec(l) + 6^0(2) + 0(63), ^ ^ ^-^2^^ ^ Q(^3)^ 

The zero-order solution is spanned by unit vectors ej at the j-th component that correspond 
to the node {n,m) G S^^^: 



E 

The first-order correction c^^^ takes the form: 



' 0, (n,m) e 

^ 0, (n,m) ^ S^USW 



where the sum defined in (|A.3|) . At the second-order in e, we find a set of non-trivial 

equations at the nodes (n, m) G S^^^: 

«i +«j+2 + 2(aj+i + aj-i) = X2CXj, j = 1,2,3,4. 

The reduced eigenvalue problem has a double zero eigenvalue and two non-zero eigenvalues 
—2 and 6. Two zero eigenvalues of the problem H5.2|l persist at all orders of e, because of the 
exact solutions: cj,„ = ^n,m and c^„j = ^n,m- ■ 

We note that the pairs of eigenvalues Ai^2) ^3^4, and Ag^io continue the eigenvalues of the vortex 
cross (1, 1) from /3 7^ 1 to /3 = 1. The pairs of eigenvalues As^g and Xj^s match with the zero 
0(e) corrections to the corresponding eigenvalues of the vortex cross (1, 1) for /? 7^ 1. Finally, 
the pair of non-zero eigenvalues Aii^i2 for /? 7^ 1 is forced to remain at the origin for /? = 1 
due to the polarization-rotation symmetry. 
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We can now specify how many purely imaginary eigenvalues A have negative Krein signature. 
When P = 1, there are four negative and twelve zero eigenvalues of W^''^ for the limiting solution 
(|2.7|1 . Out of the twelve zero eigenvalues, three small negative eigenvalues bifurcate in the 
subspace for components {an,m-,bn,m) ■, two small positive and two small negative eigenvalues 
bifurcate in the subspace for components (c^m)C^m) ^"^^ ^'^^ eigenvalues remain at zero as 
e 7^ 0. The total number of negative eigenvalues is reduced by one symmetry constraint^, 
such that eight negative eigenvalues in a constrained subspace match four pairs of imaginary 
eigenvalues with negative Krein signature. The only pair of purely imaginary eigenvalues with 
positive Krein signature is the pair As^g that is related to the two small positive eigenvalues 
in the subspace for components mi '^n m) ■ 



5.2 Eigenvalues of the (1,-1) vortex cross 



Since the stability problem ()4.1|) has no block-diagonalization for the (1,-1) vortex cross, the 
results of the second-order LS reductions give only two pairs of purely imaginary eigenvalues 
Ai_2 and A3^4. We shall study the eigenvalues of the fourth-order LS reduction by using the 
symbolic computation package based on Wolfram's Mathematica. In order to prepare for 
symbolic computations, we note that the eigenvalues of 'H(e) in the case (1,-1) are exactly 
the same as eigenvalues of 'H(e) in the case (1, 1), due to the equivalent transformation of the 
vector in the eigenvalue problem TC(€)(f = 'j(p: 



(1,-1) 







\ 






COS S 




bn,m 













C+ 











v 


r~ 


) 




v 


— sin (5 



s'md \ 

COS 5 sin 5 

— sin 5 COS 5 

COS (5 J 



As a result of this transformation, we immediately find the five- dimensional kernel of 'H(e) for 
e 7^ 0, which can be spanned as follows: 



I COS (5 $„,m ^ 

- COS 8 

- sin (5 $„^„i 
sin (5 



I -sin (5 ^ 



V cos(5$„^„ 



( \ ( -sin(5$„,„, ( 




\ cos 5 







- sin (5 
cos (5 $„,„i 




V 



- sin (5 

cos 5 ^n,m 



(5.3) 

for (n, m) E 1?. Algebraic multiplicity of zero eigenvalue for e 7^ is defined by the solution of 
the inhomogeneous equation TC{e)(p = 2ia(p, which is equivalent to the projection equations 



(n,m) 



0, 



j = 1,2,3,4,5, 



*The Hessian matrix related to two gauge symmetry constraints has a zero eigenvalue for /3 
positive eigenvalues are counted in a reduction of the negative index of 7i(e). 



1, while only 
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where </? is spanned by five eigenvectors </?j in the decomposition (|4.7jl . Solving this system 
of hnear equations, we have found under the non-degeneracy condition ()5.1() that there is a 
one-parameter solution of the inhomogeneous system for 5 7^ ^ and a three-parameter solution 
for 6 = J. Thus, the zero eigenvalue has algebraic multiplicity six for S ^ j and eight for 

In the limit e = 0, when 7i^^^ = 'H(O), we construct explicitly three sets of linearly independent 
eigenvectors of H^''^: 



/ cos 5 e'^J \ 
— cos 6 e~*^J 
sin 6 e~'^^ 
y — sin (5 e*^J J 



( cos 5 e*^J \ 



cos 5 e 



( sin b e*^J \ 



^ sin b e*^J y 



V 



sin (5 e" 
- cos b e 



cos (5 e*^J 



(5.4) 



Only the set of eigenvectors ej generates the set of generalized eigenvectors of the problem 
7^(0) ej = 2iaej, where 

/ cos 5 e*^J \ 
cos 5 e~'^^ 



sin 5 e 



\ 



sin 5 J 



(5.5) 



Thus, the zero eigenvalue of TL^^^ has algebraic multiplicity sixteen and geometric multiplicity 
twelve. Two pairs of purely imaginary eigenvalues of negative Krein signatures bifurcate at 
the second-order LS reductions as 

Ai,2,A3,4 = ±2ie + 0(e2). 

In order to study bifurcations of non-zero eigenvalues at the fourth-order LS reductions, we 
consider the extended perturbation series (|3.8() for and A with Ai = and 



j=i j=i j=i 



(5.6) 



Performing computations symbolically, we have twelve homogeneous equations at the order of 
O(e^) for twelve variables {cj,d^ ,d~), j = 1,2,3,4, which can be converted and simplified to 
the following determinant equation: 

7I + 4(1 + 4 cos 4(5)72 + 36 = 0, 



^A^. By using the inverse relation A2 = ±1/272 and finding the roots for 72 



where 72 — 2 

explicitly, we obtain four small pairs of eigenvalues with asymptotic approximations: 



^5,6 
^7,8 



±2ieVl + 4cos4(5- v^8(cos 4<5 + cos 86) + O(e^), 



±2ieVl + 4cos4(5+ v^8(cos45 + cos8(5) + 0(e^). 



17 



When (5 = or (5 = ^ , we obtain the same pairs of purely imaginary eigenvalues as in the case 
(1, 1) (see Proposition When 6 = j, we obtain two degenerate pairs of real eigenvalues 

A5,6,A7,8 = ±2^/3e2 + 0(e=^). 

The instability domain is found analytically from the condition that complex- valued roots for 
72 coalesce and become a double negative root. This happens when cos(4(5) + cos(85) = 0, 
which is solved on the interval 5 G [O, f ] ed, 6 = and 5 = Thus, the instability domain 
of the (1, —1) vortex cross in the case /3 = 1 is bounded by the interval 5 G (j^, ||). 

In order to capture the remaining pair of non-zero eigenvalues Ag^io, we shall reorder the 
perturbation series expansions and to move the last two sums in the decomposition H5.6() to the 
order of O(e^), while the coefficients of the vector c = (ci, C2, 03,04)^ should be projected to the 
vector p2 = (—1,1,-1,1) of the kernel of A^2) such that c = xip2- Performing computations 
symbolically, we have twelve homogeneous equations at the order of 0(e'^) for eight variables 
in the vectors d''' and dj and the coordinate xi. The homogeneous system is satisfied with 
the choice d+ = and = X2P21 where X2 is another coordinate. The coordinates (xi,X2) 
solve a homogeneous system with the determinant equation A| = — 16cos^(25). Therefore, a 
small pair of purely imaginary eigenvalues of negative Krein signatures has the asymptotic 
approximation: 

A9,io = ±4^2 cos(25) + O(e^). 

When (5 = and 6 = ^, the pair Ag^io matches to that in the case (1,1) (see Proposition 
15. If) . When S = j, the pair remains at the origin as it follows from the study of algebraic 
multiplicity of zero eigenvalue. According to the count of negative eigenvalues, the total 
number of negative eigenvalues of 'H(e) for small e reduced by one symmetry constraint is 
eight. These eigenvalues match two pairs of imaginary eigenvalues Ai^2 and A3^4 and two real 
positive eigenvalues As^g and Xj^s^. 

Asymptotic and numerical approximations of small eigenvalues A for small values of e for 
uj = P = 1 and S = J are shown on Figure |SJ The left plot corresponds to the vortex 
pair (1,1), while the right plot corresponds to the vortex pair (1,-1). We can see that the 
(1,1) vortex cross is linearly stable in the anti-continuum limit, according to the results of 
Proposition 15.11 On the other hand, the (1,-1) vortex cross become unstable because of 
the a double pairs of real eigenvalues As^g = -^7,8- The other double pair of purely imaginary 
eigenvalues remains double for all e > 0, such that Ai^2 = A3,4. Therefore, the stability changes 

^Another solution exists at 5 = ^ but it corresponds to the case when complex-valued roots coalesce and 
become a double positive root for 72. 

^Eigenvalues As.e and Ay.g are real only in the case S = j. For other values of S, these eigenvalues are either 
complex- valued or purely imaginary. The count is not affected, since two real eigenvalues are equivalent to four 
complex eigenvalues which may coalesce due to the inverse Hamilton-Hopf bifurcation to two pairs of purely 
imaginary eigenvalues with positive and negative Krein signatures. 
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Figure 5: Eigenvalues of the vector vortex cross with to = (3 = 1 and S = j versus e. Left: 
(1,1). Right: (1,-1). The sohd hnes show the full numerical results, while the dashed lines 
show the asymptotic approximations. Bolded curves show double eigenvalues. 

drastically in the case of the discrete Manakov system (that is the coupled DNLS system for 
/? = 1): the (1, 1) vortex cross is stable near the anti-continuum limit while the (1,-1) vortex 
cross is linearly unstable. 

6 Conclusions 

We have examined analytically and numerically the existence and stability of vortex cross 
configurations in the single-component and two-component DNLS equations. We have used the 
Lyapunov-Schmidt theory, to obtain the bifurcation functions and the solvability conditions 
that allow persistence of such configurations near the anti-continuum limit. Additionally, 
the theory gives analytical expressions for eigenvalues of the linearized stability problem as 
functions of the system parameters (namely, the coupling between adjacent lattice nodes e > 
and the coupling between the two components /3 > 0). 

One of the interesting recent experimental developments in the setting of BECs concerns the 
experimental and theoretical studies of spin-1 (or 3-component) states This feature, 

along with similar possibilities that could be realized in optical settings, render desirable a 
general theory for interactions of multiple components. Such studies are currently in progress. 
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A Continuation of the single-component vortex cross 



We apply the algorithm of Lyapunov-Schmidt (LS) reductions (see [5] for details) and compute 
the first few terms of the perturbation series expansions: 



(k) 
m' 



(A.l) 



k=0 



The zero-order solution ^n)m is given by H3.2|) . The first-order correction is obtained in the 
explicit form: 

0, (n,m)G5(0) 

Ep^e^^S (n,m)G5W (A.2) 
0, (n,m) ^ 5(0)U5« 

where 5^^^ is the set of adjacent nodes to the set S'^'^^ and X^P^ 6*^' is a schematic notation for 
the following solution: 



$(1) 



(1) 
I 



f,tOi gi02 + + f,t04^ (^n, m) = (0, 0) 
[ e*^^ (n,m) ={(-2,0); (0,-2); (2,0); (0,2)} 



(A.3) 



The index j enumerates nodes in the set S^'^^ that are adjacent to the nodes in the set S^^^ 
listed in the figured brackets of ()A.3|) . No non-trivial bifurcation equations arise at the first- 
order reductions, i.e. the first-order correction to the bifurcation function g*^^)(0) is zero, 
where 6 = {61,62,03,04^) and notations of |9| are used. The second-order correction is found 
in the form: 



$(2) 



0, 

E 



n, m 



(2) AQi 



G5« 



n, ra 



(A.4) 



where 



-2s 



(2) 



4 + 2 cos(6'j+i - 6.j) + 2 cos(6'j_i - 6.j) + cos(0j+2 - 6j). 

The set ^'^^-^ contains outward adjacent nodes to the set 5*^^-'\{(0, 0)} and Y^i^ ^^^^ ^ 
schematic notation for the following solution: 



(2) 



r 2e*^^ +e^^^+i, (n,m) ={(-2,-1); (1,-2); (2,1); (-1,2)} 
!^^' = | e^e, +2ee,+i^ (n, m) = {(-1, -2); (2, -1); (1, 2); (-2, 1)} 
' [ e^^^ (n,m) = {(-3,0); (0,-3); (3,0); (0,3)} 

The second-order corrections to the bifurcation function take the form: 

gf = 2 sin(0j - 6j+i) + 2 sin(0, - 0,_i) + sin(0,- - 6^+2), j = 1, 2, 3, 4. 



(A.5) 



(A.6) 
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The bifurcation equations g^'^\0) = are satisfied with the one-parameter family of asym- 
metric vortices: 

9i = 0, 02 = 6, 03 = 7r, 6^ = 7: + e, (A.7) 

where 6 G (0, vr). When = ^, the family HA.7|) reduces to the vortex cross configuration 
(|2.12|) . The Jacobian matrix M.2 of the second-order bifurcation function g(^^(0) is obtained 
by differentiation of g(^) in 6. At the family of asymmetric vortices (|A.7() . the Jacobian matrix 
A^2 takes the form: 

/ -1 -2 cos 6* 1 2 cos 6* \ 
-2 cos 6* -1 2 cos 6* 1 
1 2 cos 61 -1 -2 cos 61 

\^ 2 cos 6* 1 -2 cos 6* -1 j 



It has two zero eigenvalues and two non-zero eigenvalues —2 it 4 cos 9. In the case of the vortex 
cross {6 = it has two zero eigenvalues and two negative eigenvalues —2. The third-order 
correction satisfies the inhomogeneous equation, 



fl_2|$(0) |2)$(3) _ $(0)2^(3) ^^(2) ^(2) ^(2) ,^(2) ,|$(1)|2^(1 

^\^n,m\ J^n,m ^n,'m^n,m n+l,m ' n—l,m ' n,m+l ' n,m—l ' I n, ml n, 



,(2) 



,(2) 



1) 



where we have shorten nonlinear terms, 

The third-order correction is found in the form: 



for all (n, m) G 



$(3) 



0, (n,m)G5(°) 
0, (n,m)G5(2) 
^ 0, (n, m) i 5(0) u 5(1) u 5(2) u 5(3) 



(A. 



where the sum X^P'* sp^e*^' is defined similarly to the sum (TOl) . the sum is '^ot used 

for further computations, and the sum is defined as follows: 



(1,3) 



5e«ei +e*^^+i +e*^^-i, (n,m) = {(-2, 0); (0, -2); (2, 0); (0, 2)} 



(A.9) 



No non-trivial bifurcation equations arise at the third-order reductions, i.e. g('^^(0) = 0. The 
fourth-order correction satisfies the inhomogeneous equation, 

fl_2l$(°) Pl*^'^^ _$(0)2|,(4) ^2l$(^)P$(°^ +$(2)2$(0) i$(3), +$(3) ^(3) ,^(3) 

V I n. ml / n,m n,m n,m I n, ml n,m * n,m n,m ' Ti+l,m ' 7i — l,m ' n,m+l ' n,m — 1" 

Solving the inhomogeneous equation for the third-order corrections, we obtain the bifurcation 
equations at the fourth order of LS reductions in the form: 

gf^ = (4 + 2 cos(6ij+2 - Oj+i) + 2 cos(6ij - Oj+i) + cos(e'j_i - 6ij+i)) sin(6'j+i - 61^) 
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+ (4 + 2 cos(6'j_2 - 6*^-1) + 2 cos(6'j - 6j^i) + cos(6'j+i - 
+ ^ (4 + 2 cos(6lj_i - 6lj_2) + 2 cos(6lj+i - 61^+2) + cos(0j 

+ i(4 + 2cos((.,_ 
+ 2(1 + cos(6lj+i 

+ 2(2 + C0S(6I2 - 6*1) + COS(03 - Oi) + 008(6*4 - 6*1) + cos(6i3 - O2) + C0S(6'4 - 02) + cos{94 - 6I3)) 
X (sin(6'j - Oj+i) + sin(6lj - 6lj_i) + sin(6'j - 6'j+2)) 



',_i))sin(0,_i-0,) 
- 0j+2))sin(6lj+2 - 6lj) 

(4 + 2cos(6lj+i - 0J+2) + 2cos(6lj_i - 6^.2) + cos(6lj - 0j_2)) sin(6'j_2 - ^j) 
j)) sin(0j - Oj+i) + 2(1 + cos(6lj_i - 0^)) 8111(6*^ - 9j^i) 



4sin(0j - 



For the asymmetric vortex, we have 



5f) = (-l)^2sin(2e), 



1,2,3,4. 



The Jacobian matrix has two zero eigenvalues with orthogonal eigenvectors: 

1 1 

1 ' -1 



Pi 



V 1 J 

It is clear that the vector g^^^^ = 2 sin(20)p2 is not orthogonal to the eigenvector p2 of the kernel 
of Ai2, unless 9 = {0, f ,vr}. By Proposition 2.10 in j^], the family of asymmetric vortices 
(|A.7|) terminates at the fourth-order reduction. The exceptional cases include discrete solitons 
for 9 = {0, tt} and the vortex cross at ^ = ^. In order to consider persistence of the vortex 
cross, we compute the Jacobian matrices Ai2 and 7W4 from the bifurcation functions g*-^-* and 
g^^) explicitly. 



/ -1 

1 

V 





1 



1 \ 

1 

-1 

-I J 



( 3 
2 



-7 



V 



2 

3 2 
2 3 
-7 2 



7 



2 \ 

-7 
2 

3 / 



Since A^4Pi = and A44P2 / 0, the zero eigenvalue of M2 with the associated eigenvector p2 
bifurcates. By Proposition 2.9 in 9 , this implies that the family of the vortex cross is continued 
from the anti-continuum limit uniquely up to the rotational transformation 0^0 + ^oPi that 
corresponds to the gauge symmetry of the dNLS equation ()3.1|) . Proposition 13.11 is hence 
proved. 

Small eigenvalues of the linearized Jacobian matrix 'H(e) are defined by an extended eigenvalue 
problem for the Jacobian matrices M2 and A^4, 



{e^M2 + e^MA + 0{e^))c = ^c. 
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There exist four eigenvalues of the extended problem which admit the asymptotic approxima- 
tions, 

7i,2 = -2e2 + 0(e4), 73 = -Se^ + 0(e6), 74 = 0. 

The eigenvalue 73 is obtained by the perturbation theory for the zero eigenvalue of A^2 asso- 
ciated with the eigenvector p2 (orthogonal to the eigenvector pi): 

hm e 73 = — — = -8. 

(P2,P2) 

References 

[1] J.W. Fleischer, G. Bartal, O. Cohen, T. Schwartz, O. Manela, B. Preedman, M. Segev, 
H. Buljan, and N.K. Efremidis, "Spatial photonics in nonlinear waveguide arrays". Opt. 
Express 13, 1780 (2005). 

[2] V.A. Brazhnyi and V.V. Konotop, "Theory of nonlinear matter waves in optical lattices". 
Mod. Phys. Lett. B 18, 627 (2004) 

[3] P.G. Kevrekidis and D.J. Frantzeskakis, "Pattern forming dynamical instabilities in Bose- 
Einstein condensates", Mod. Phys. Lett. B 18, 173 (2004) 

[4] O. Morsch and E. Arimondo, in Dynamics and Thermodynamics of Systems with Long- 
Range Interactions, T. Dauxois, S. Ruffo, E. Arimondo and M. Wilkens (Eds.), Springer 
(Berlin 2002), pp. 312-331. 

[5] P.G. Kevrekidis, K.0. Rasmussen and A.R. Bishop, "The discrete nonlinear Schrodinger 
equation: a survey of recent results". Int. J. Mod. Phys. B 15, 2833 (2001) 

[6] G.L. Alfimov, P.G. Kevrekidis, V.V. Konotop and M. Salerno, "Wannier functions anal- 
ysis of the nonlinear Schrodinger equation with a periodic potential", Phys. Rev. E 66, 
046608 (2002) 

[7] Yu.S. Kivshar and M. Peyrard, Phys. Rev. A 46, 3198 (1992). 

[8] D.E. Pelinovsky, P.G. Kevrekidis, and D.J. Frantzeskakis, "Stability of discrete solitons 
in nonlinear Schrodinger lattices", nlin.PS/0410005[ Physica D, in press (2005). 

[9] D.E. Pelinovsky, P.G. Kevrekidis, and D.J. Frantzeskakis, "Persistence and stability of 
discrete vortices in nonlinear Schrodinger lattices" , |nlin. PS/041 1016[ Physica D, in press 
(2005). 

[10] J.W. Fleischer, T. Carmon, M. Segev, N.K. Efremidis and D.N. Christodoulides, "Obser- 
vation of discrete solitons in optically induced real time waveguide arrays", Phys. Rev. 
Lett. 90 023902 (2003) 



23 



[11] H. Martin, E.D. Eugcnicva, Z. Chen and D.N. Christodoulides, "Discrete solitons and 
soliton-induced dislocations in partially coherent photonic lattices", Phys. Rev. Lett. 92 
123902 (2004) 

[12] J. Yang, I. Makasyuk, A. Bezryadina and Z. Chen, "Dipole solitons in optically-induced 
two-dimensional photonic lattices" , Opt. Lett. 29, 1662 (2004) 

[13] Z. Chen, H. Martin, E.D. Eugenieva, J. Xu and A. Bezryadina, "Anisotropic enhancement 
of discrete diffraction and formation of two-dimensional discrete-soliton trains", Phys. 
Rev. Lett. 92 143902 (2004) 

[14] J. Yang, I. Makasyuk, P.G. Kevrekidis, H. Martin, B.A. Malomed, D.J. Frantzeskakis, 
and Z. Chen, " Necklace- like solitons in optically induced photonic lattices", Phys. Rev. 
Lett. 94, 113902 (2005) 

[15] Z. Chen, I. Makasyuk, A. Bezryadina and J. Yang, "Observation of two-dimensional 
lattice vector solitons". Opt. Lett. 29 1656 (2004) 

[16] D.N. Neshev, T.J. Alexander, E.A. Ostrovskaya, Yu.S. Kivshar, H. Martin, L Makasyuk 
and Z. Chen, "Observation of discrete vortex solitons in optically induced photonic lat- 
tices", Phys. Rev. Lett. 92, 123903 (2004) 

[17] J.W. Fleischer, G. Bartal, O. Cohen, O. Manela, M. Segev, J. Hudock and D.N. 
Christodoulides, "Observation of vortex-ring "discrete" solitons in 2D photonic lattices", 
Phys. Rev. Lett. 92 (2004) 123904 

[18] B.A. Malomed and P.G. Kevrekidis, "Discrete vortex solitons", Phys. Rev. E 64, 026601 
(2001) 

[19] J. Yang and Z. Musslimani, "Fundamental and vortex solitons in a two-dimensional op- 
tical lattice". Opt. Lett. 23, 2094 (2003) 

[20] P.G. Kevrekidis, B.A. Malomed, Z. Chen and D.J. Frantzeskakis, "Stable higher-order 
vortices and quasi- vortices in the discrete nonlinear Schrodinger equation" , Phys. Rev. E 
70, 056612 (2004). 

[21] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Sengstock, A. Sanpera, G.V. Shlyapnikov 
and M. Lewenstein, "Dark solitons in Bose-Einstein condensates", Phys. Rev. Lett. 83, 
5198 (1999) 

[22] J. Denschlag, J.E. Simsarian, D.L. Feder, C.W. Clark, L.A. Collins, J. Cubizolles, L. 
Deng, E.W. Hagley, K. Helmerson, W.P. Reinhardt, S.L. Rolston, B.L Schneider and 
W.D. Phillips, "Generating solitons by phase engineering of a Bose-Einstein condensate". 
Science 287, 97 (2000) 



24 



[23] B.P. Anderson, P.C. Haljan, C.A. Regal, D.L. Fcdcr, L.A. Collins, C.W. Clark and E.A. 
Cornell, "Watching dark solitons decay into vortex rings in a Bose-Einstein condensate", 
Phys. Rev. Lett. 86, 2926 (2001) 

[24] K. E. Strecker, G.B. Partridge, A.G. Truscott and R.G. Hulet, "Formation and propaga- 
tion of matter-wave soliton trains". Nature 417, 150 (2002) 

[25] L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. CubizoUes, L.D. Carr, Y. Castin 
and C. Salomon, "Formation of a matter- wave bright soliton". Science 296, 1290 (2002) 

[26] B. Eiermann, Th. Anker, M. Albiez, M. Taglieber, P. Treutlein, K.-P. Marzhn and M.K. 
Oberthaler, "Bright Bose-Einstein gap solitons of atoms with repulsive interaction" , Phys. 
Rev. Lett. 92, 230401 (2004). 

[27] M. Grcincr, I. Bloch, O. Mandcl, T.W. Hansch and T. Esslingcr, "Bosc-Einstcin Conden- 
sates in ID and 2D optical lattices", Appl. Phys. B 47, 769 (2001) and "Exploring phase 
coherence in a 2D lattice of Bose-Einstein condensates", Phys. Rev. Lett. 87, 160405 
(2001). 

[28] J. Meier, J. Hudock, D. Christodoulides, G. Stegeman, Y. Silberberg, R. Morandotte and 
J.S. Aitchison, "Discrete vector solitons in Kerr nonlinear waveguide arrays", Phys. Rev. 
Lett. 91, 143907 (2003) 

[29] J. Meier, J. Hudock, D.N. Christodoulides, G.L Stegeman, H.Y. Yang, G. Salamo, R. 
Morandotti, J.S. Aitchison and Y. Silberberg, "Excitation of strongly confined scalar and 
vector self-trapped beams in one-dimensional arrays of Kerr-nonlinear channel waveg- 
uides", J. Opt. Soc. Am. B 22, 1432 (2005). 

[30] C. J. Myatt, E.A. Burt, R.W. Christ, E.A. Cornell and C.E. Wieman, "Production of 
two overlapping Bose-Einstein condensates by sympathetic cooling", Phys. Rev. Lett. 78, 
586 (1997) 

[31] D.S. Hall, M.R. Matthews, J.R. Enshcr. C.E. Wicman and E.A. Cornell, "Dynamics of 
component separation in a binary mixture of Bose-Einstein condensates", Phys. Rev. 
Lett. 81, 1539 (1998) 

[32] D.M. Stamper-Kurn, M.R. Andrews, A. P. Chikkatur, S. Inouye, H.-J. Miesner, J. Stenger 
and W. Ketterle, "Optical confinement of a Bose-Einstein condensate", Phys. Rev. Lett. 
80, 2027 (1998) 

[33] G. Modugno, G. Ferrari, G. Roati, R.J. Brecha, A. Simoni and M. Inguscio, "Bose- 
Einstein condensation of potassium atoms by sympathetic cooling". Science 294, 1320 
(2001) 



25 



[34] M. Mudrich, S. Kraft, K. Singer, R. Grimm, A. Mosk and M. Weidemiiller, "Sympathetic 
cooling with two atomic species in an optical trap", Phys. Rev. Lett. 88, 253001 (2002). 

[35] F. Lederer, S. Darmanyan and A. Kobyakov in Spatial Optical Solitons, edited by S. Trillo 
and W.E. Torruellas (Springer- Verlag, New York, 2001) 

[36] M.J. Ablowitz and Z.H. Musslimani, "Discrete vector spatial solitons in a nonlinear waveg- 
uide array", Phys. Rev. E 65, 056618 (2002) 

[37] J. Hudock, P. G. Kevrekidis, B. A. Malomed, and D. N. Christodoulides, "Discrete vector 
solitons in two-dimensional nonlinear waveguide arrays: solutions, stability and dynam- 
ics", Phys. Rev. E 67, 056618 (2003) 

[38] R.A. Vicencio, M.I. Molina and Yu.S. Kivshar, "All-optical switching and amplification 
of discrete vector solitons in nonlinear cubic birefringent waveguide arrays". Opt. Lett. 
29, 2905 (2004) and "Polarization instability, steering and switching of discrete vector 
solitons", Phys. Rev. E 71, 056613 (2005) 

[39] A.S. Desyatnikov, D. Michalache, D. Mazilu, B.A. Malomed, C. Denz and F. Lederer, 
"Two-dimensional solitons with hidden and explicit vorticity in bimodal cubic-quintic 
media", Phys. Rev. E 71, 026615 (2005) 

[40] F. Ye, J. Wang, L. Dong and Y.-P. Li, "Suppression of modulational instability of ring 
vector solitons". Opt. Comm. 230, 219 (2004) 

[41] D.E. Pelinovsky and Yu.S. Kivshar, "Stability criterion for multicomponent solitary 
waves", Phys. Rev. E 62, 8668-8676 (2000). 

[42] J.M. Higbie, L.E. Sadler, S. Inouye, A.P. Chikkatur, S.R. Leslie, K.L. Moore, V. Savalh 
and D.M. Stamper-Kurn, "Direct nondestructive imaging of magnetization in a spin-1 
Bose-Einstein gas", Phys. Rev. Lett. 95, 050401 

[43] H. Saito and M. Ueda, "Spontaneous magnetization and structure formation in a spin-1 
ferromagnetic Bose-Einstein condensate", cond-mat /0504398j 



26 



